Proteomics Report: DIA Benchmark Experiment

DE Analysis

Author

Andrea

Experiment

Yeast as background and USP2 proteins spiked in different concentration. This Experiment is designed for DIA benchmarking of different workflow using DIA-NN.

Import data

Reading DIA data from data/report.tsv using Precursor.Translated intensities from DIA-NN.

All the precursors with Q-value < 0.01 are used for the analysis.

Creating Q-Features object and run pre processing

The pre- processing steps are:

  • Filtering precursor with more than 3 samples with no missing intensities across all the samples
  • Select only proteotypic peptides
  • Select proteins with at least 3 peptides

Data Missing Analysis

The plot shows the completness of the experiments at precursor and summarized level. On the the x-axis peptides/proteins are ordered based on the increasing NA content.

Normalization

Raw intensities from DIA-NN are processed as follow:

  • Log 2 transformation
  • Normalized across all the samples using quantiles method
  • Summarization at protein level using robustSummary function

PCA

DE Analysis

Using MSqRob2.

  GroupA   GroupB   GroupC   GroupD   GroupE 
15.13545 14.78455 14.88787 15.01882 15.01860 

QC plots

Group Comparison

This volcano plot summarizes the differential expression landscape in the comparison between the two groups

Differential Expressed Proteins 

These bar plots summarizes the number of significantly upregulated/downregulated number of miRNAs based on different adjusted p-values (selected adjusted p-values are 0.001, 0.01, 0.05, and 0.1 - see facet headers) and log2 fold-change thresholds (on the x-axis) used to define the significance levels.

This volcano plot summarizes the differential expression landscape in the comparison between the two groups

Differential Expressed Proteins 

These bar plots summarizes the number of significantly upregulated/downregulated number of miRNAs based on different adjusted p-values (selected adjusted p-values are 0.001, 0.01, 0.05, and 0.1 - see facet headers) and log2 fold-change thresholds (on the x-axis) used to define the significance levels.

Session Info

R version 4.3.3 (2024-02-29 ucrt)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 19045)

Matrix products: default


locale:
[1] LC_COLLATE=English_Belgium.utf8  LC_CTYPE=English_Belgium.utf8   
[3] LC_MONETARY=English_Belgium.utf8 LC_NUMERIC=C                    
[5] LC_TIME=English_Belgium.utf8    

time zone: Europe/Brussels
tzcode source: internal

attached base packages:
[1] stats4    stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
 [1] reshape2_1.4.4              pheatmap_1.0.12            
 [3] biomaRt_2.58.2              gridExtra_2.3              
 [5] ggVennDiagram_1.5.2         hrbrthemes_0.8.7           
 [7] factoextra_1.0.7            MSnbase_2.28.1             
 [9] ProtGenerics_1.34.0         mzR_2.36.0                 
[11] Rcpp_1.0.12                 DT_0.33                    
[13] plotly_4.10.4               here_1.0.1                 
[15] magrittr_2.0.3              lubridate_1.9.3            
[17] forcats_1.0.0               stringr_1.5.1              
[19] purrr_1.0.2                 readr_2.1.5                
[21] tidyr_1.3.1                 tibble_3.2.1               
[23] ggplot2_3.5.0               tidyverse_2.0.0            
[25] dplyr_1.1.4                 visdat_0.6.0               
[27] readxl_1.4.3                msqrob2_1.10.0             
[29] QFeatures_1.12.0            MultiAssayExperiment_1.28.0
[31] SummarizedExperiment_1.32.0 Biobase_2.62.0             
[33] GenomicRanges_1.54.1        GenomeInfoDb_1.38.8        
[35] IRanges_2.36.0              S4Vectors_0.40.2           
[37] BiocGenerics_0.48.1         MatrixGenerics_1.14.0      
[39] matrixStats_1.3.0          

loaded via a namespace (and not attached):
  [1] splines_4.3.3           later_1.3.2             bitops_1.0-7           
  [4] filelock_1.0.3          cellranger_1.1.0        preprocessCore_1.64.0  
  [7] XML_3.99-0.16.1         lifecycle_1.0.4         rstatix_0.7.2          
 [10] doParallel_1.0.17       rprojroot_2.0.4         lattice_0.22-6         
 [13] MASS_7.3-60.0.1         crosstalk_1.2.1         backports_1.4.1        
 [16] sass_0.4.9              limma_3.58.1            rmarkdown_2.26         
 [19] jquerylib_0.1.4         yaml_2.3.8              httpuv_1.6.15          
 [22] RColorBrewer_1.1-3      MsCoreUtils_1.14.1      DBI_1.2.2              
 [25] minqa_1.2.6             abind_1.4-5             zlibbioc_1.48.2        
 [28] AnnotationFilter_1.26.0 RCurl_1.98-1.14         rappdirs_0.3.3         
 [31] gdtools_0.3.7           GenomeInfoDbData_1.2.11 ggrepel_0.9.5          
 [34] crul_1.4.2              ncdf4_1.22              codetools_0.2-20       
 [37] DelayedArray_0.28.0     xml2_1.3.6              tidyselect_1.2.1       
 [40] farver_2.1.1            httpcode_0.3.0          lme4_1.1-35.3          
 [43] BiocFileCache_2.10.2    jsonlite_1.8.8          iterators_1.0.14       
 [46] systemfonts_1.0.6       foreach_1.5.2           tools_4.3.3            
 [49] progress_1.2.3          glue_1.7.0              BiocBaseUtils_1.4.0    
 [52] Rttf2pt1_1.3.12         SparseArray_1.2.4       xfun_0.43              
 [55] withr_3.0.0             BiocManager_1.30.22     fastmap_1.1.1          
 [58] boot_1.3-30             fansi_1.0.6             digest_0.6.35          
 [61] timechange_0.3.0        R6_2.5.1                mime_0.12              
 [64] colorspace_2.1-0        RSQLite_2.3.6           utf8_1.2.4             
 [67] generics_0.1.3          fontLiberation_0.1.0    data.table_1.15.4      
 [70] prettyunits_1.2.0       httr_1.4.7              htmlwidgets_1.6.4      
 [73] S4Arrays_1.2.1          pkgconfig_2.0.3         gtable_0.3.4           
 [76] blob_1.2.4              impute_1.76.0           XVector_0.42.0         
 [79] htmltools_0.5.8.1       carData_3.0-5           fontBitstreamVera_0.1.1
 [82] MALDIquant_1.22.2       clue_0.3-65             scales_1.3.0           
 [85] png_0.1-8               knitr_1.46              tzdb_0.4.0             
 [88] nlme_3.1-164            curl_5.2.1              nloptr_2.0.3           
 [91] cachem_1.0.8            parallel_4.3.3          extrafont_0.19         
 [94] AnnotationDbi_1.64.1    mzID_1.40.0             vsn_3.70.0             
 [97] pillar_1.9.0            grid_4.3.3              vctrs_0.6.5            
[100] pcaMethods_1.94.0       ggpubr_0.6.0            promises_1.3.0         
[103] car_3.1-2               dbplyr_2.5.0            xtable_1.8-4           
[106] cluster_2.1.6           extrafontdb_1.0         evaluate_0.23          
[109] cli_3.6.2               compiler_4.3.3          rlang_1.1.3            
[112] crayon_1.5.2            ggsignif_0.6.4          labeling_0.4.3         
[115] affy_1.80.0             plyr_1.8.9              stringi_1.8.3          
[118] viridisLite_0.4.2       BiocParallel_1.36.0     munsell_0.5.1          
[121] Biostrings_2.70.3       lazyeval_0.2.2          fontquiver_0.2.1       
[124] Matrix_1.6-5            hms_1.1.3               bit64_4.0.5            
[127] gfonts_0.2.0            KEGGREST_1.42.0         statmod_1.5.0          
[130] shiny_1.8.1.1           broom_1.0.5             igraph_2.0.3           
[133] memoise_2.0.1           bslib_0.7.0             affyio_1.72.0          
[136] bit_4.0.5